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APPUCATION  OF  RAINFALL-RUNOFF 
SIMULATION  FOR  FLOOD  FORECASTlNGl 

John  Peters^ 


INTBGDainON 

Hydrologic  simulation  models  provide  a  means  for  extending  the  lead  time 
associated  with  flood  warnings.  Following  a  brief  discussion  of  warning  ob¬ 
jectives  and  approadies  to  short-term  hydrologic  forecasting,  characteristics 
of  rainfall-runoff  models  for  forecast  ai^llcatlons  are  discussed.  HEC-1  and 
HEC-2  are  described  as  Illustrations  of  models  that  can  be  applied  to  develop 
warning  criteria;  HECIF  and  the  Sacramento  Method  are  described  as  Illustra¬ 
tions  of  models  for  real-time  application.  Finally,  aspects  of  model  selec- 
tlc»i  and  use  are  discussed. 


HABHINS  OeJBCnVBS 

The  value  of  a  flood-threat  recognition  system  depends  to  a  large  extent 
on  the  lead  time  it  provides  for  Issuing  wartungs,  aiabllng  evacuation,  etc. . 

A  minimum  lead  time  must  be  provided  for  a  system  to  be  practically  useful. 

The  lead  time  that  Is  potentially  achievable  depends  on  (1)  the  spatial  and 
temporal  characteristics  of  storm  rainfall  and  the  ability  to  sample/forecast 
these,  (2)  rainfall-runoff  response  characteristics  of  the  watershed  and  the 
ability  to  simulate  these,  and  (3)  the  time  required  to  recognize  and  evalu¬ 
ate  the  flood  threat  and  take  appropriate  acticxi.  The  value  of  a  warning 
system  depends  also,  of  course,  on  its  reliability.  Consider  Figure  1,  vdiich 
illustrates  aspects  of  reliability.  Sets  of  storm  events  are  labeled  {A}, 

{B},  {C}  and  {D},  »diere: 

{A}  =  storm  events  thal  cause  flooding 

{B}  =  storm  events  that  do  not  cause  flooding 

{C}  =  storm  events  that  cause  flooding  but  for 

which  warnings  eure  not  issued 

{D}  =  storm  events  that  do  not  cause  flooding  but 

for  warnings  are  issued 


1  Paper  presented  at  the  United  States  -  Republic  of  China  Workshop  on  Natural 
Disaster  Reduction,  Taipei,  Taiwan,  June  24-26,  1993. 

2  Senior  Engineer,  US  Army  Corps  of  Engineers,  Hydrologic  Engineering  Center, 
Davis,  California,  95616,  USA. 


Fi£ure  1.  BSLIABIUTf  OF  ELOGD  HAEMDES 


The  goal  of  a  warning  system  is  to  minimize  both  {C}  and  {D}.  Events 
from  {C}  can  cause  damage  and  loss  of  life  that  could  possibly  be  prevented; 
events  from  {D}  increase  the  likelihood  that  future  warnings  will  be  ignored. 
Alternative  warning  systems  will  be  reflected  by  different  configurations  of 
{C}  and  {D}. 

Ihe  basis  for  a  warning  can  range  frcxa  measured  river  stage  (elevation) 
at  an  index  gage  to  results  of  a  rainfall-runoff  simulation  that  Incorporates 
recent  rain  data  and  estimates  of  future  rainf^dl.  Althou^  the  more  soj^is- 
ticated  warning  systems  may  provide  longer  lead  tixaes,  their  reliability  is 
not  necessarily  greater  th^  that  associated  with  simpler  systems.  Both  lead 
time  and  reliability  should  be  evaluated  vdien  analyzing  alternative  warning 
syst^Ds. 

The  tradeoff  between  lead  time  (warning  time)  and  warning  reliability 
can  be  illustrated  by  considering  a  simple  threshold-stage  method  of  warning, 
as  shown  in  Figure  2  .  The  warning  stage  is  sensed  at  location  A.  Ihe  pri¬ 
mary  flood  threat  is  downstream  at  location  B.  The  problem  is  to  choose  a 
threshold  (index)  stage  for  locatlcxi  A  such  that  vdien  that  stage  is  exceeded, 
a  warning  for  flooding  at  location  B  is  to  be  issued.  It  is  desired  that  the 
lead  time  to  prepare  for  the  flood  threat  be  as  l<»ig  as  possible.  The  lower 
the  index  stage  at  A,  presumably  the  more  lead  time  will  be  provided.  How¬ 
ever,  if  the  threshold  stage  is  too  low,  there  will  be  too  many  false  warn¬ 
ings,  so  that  genuine  warnings  will  not  be  heeded.  In  terms  of  Figure  1,  as 
{C}  becomes  smaller,  {D}  becomes  larger. 
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Figure  2.  LEAD  TIME  VS.  HftBMIMS  HKT.TABn.nT 


The  graite  in  the  lower  portion  of  Figure  2  represent  relationships  that 
could  be  developed  by  analyzing  a  set  of  historical  storm  events  (Dotson  and 
Peters,  1990).  Lead  time  is  a  variable  that  depends  on  event-specific  storm 
and  runoff  characteristics.  Storm  and  streamflow  data  from  historical  events 
provide  useful  information  for  assessing  the  magnitude  and  variability  of  lead 
time. 


AFFRQACHBS  TO  SHOBT-TQH  HYDEKXiOGIC  FDRBCASTINS 

"Short-term"  here  refers  to  forecasts  with  lead  times  of  hours  to  sev¬ 
eral  days,  as  required  for  flood  warning  purposes.  By  contrast,  long-term 
forecasts,  which  provide  lead  times  up  to  a  year  or  more,  are  useful  for  water 
management  decisions. 
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Hydrologic  models  for  short- term  forecasting  may  anploy  channel  routing, 
rainfall-runoff  simulation,  or  both.  Choice  of  a  model  type  depends  on  re¬ 
quired  forecast  lead  time,  the  response  c^iaracteristics  of  the  basin,  and  the 
scale  of  meteorological  events.  Lettenmaier  and  Wood  (Maidmait,  1993)  list 
four  cases,  as  follows: 

HawA  1  Required  lead  time  is  larger  than  the  hydrologic  response  time 
(the  sum  of  the  time  of  concentration  for  the  basin  and  the  time  of 
travel  throu^  the  river  system).  In  this  case,  a  forecast  of  precipi¬ 
tation  is  required,  because  future  precipitation  will  reach  the  forecast 
point  within  the  lead  time.  Rainfall-runoff  modeling  is  required. 

Hafw  9  Required  lead  time  is  smaller  than  the  hydrologic  response 
time,  and  the  time  of  ooncentratlcxi  is  substantially  smaller  than  the 
time  of  travel  throu^  the  river  system.  This  is  the  case  for  large 
river  systems  for  fdiich  forecasts  can  be  based  on  channel  routing  of  up¬ 
stream  observed  (gauged)  flows. 

Case  3.  Required  lead  time  is  smaller  than  the  hydrologic  response 
time,  and  the  time  of  conoentratl<»i  is  substantially  larger  than  the 
time  of  travel  throu^  the  river  system.  A  rainfall-runoff  model  is  re¬ 
quired,  but  forecasts  of  future  precipitation  are  not  required. 

Case  4.  This  case  is  cne  in  »Aiich  the  scale  of  the  meteorological  event 
is  significantly  smaller  than  the  scale  of  the  basin.  This  would  occur, 
for  example,  on  large  basins  subject  to  convective  storms.  In  this 
situation,  it  is  necessary  to  subdivide  the  basin  or  employ  a  model  that 
permits  specification  of  precipitation  as  a  distributed  input.  Also  it 
is  desirable  to  telaneter  data  from  stream  gages  on  major  tributaries. 
Rainfall-runoff  modeling  and  channel  routing  of  observed  and/or  fore¬ 
casted  tributary  flows  are  required. 

For  quick-responding  watersheds  (i.e.  Case  1  above),  lead  times  are  very 
short,  and  the  time  available  for  processing  forecasts  is  extremely  limited. 

In  such  situations,  real-time  modeling  may  not  be  practical,  and  it  may  be 
more  appropriate  to  apply  pre-established  warning  criteria  with  real-time 
observed  and  forecasted  rainfall  depths/durations  as  inputs.  Historical 
events  can  be  analyzed  to  develop  such  criteria.  The  hydrologic  models  used 
for  this  purpose  do  not  require  the  special  functionality  associated  with 
real-time  applications. 

Another  approach  for  quick-responding  watersheds  is  to  automatically 
compute  forecasts  at  frequent  intervals  (e.g. ,  every  10  minutes),  and  to  in¬ 
clude  in  those  forecasts  a  set  of  pre-specified  (fixed)  future  rainfall 
amounts,  as  illustrated  in  Figure  3.  The  amount  of  forecasted  rainfall  can 
then  be  used  to  interpolate  a  discharge  hydrograi^  (and  associated  inundated 
cuea)  from  the  most  recent  pre-coraputed  set  of  forecasts.  Automated  forecast¬ 
ing  requires  a  well  calibrated,  robust  model.  The  state-of-the-art  of  model¬ 
ing  is  such  that  generally  an  experienced  modeler  must  be  involved  in  model 
applications  and  interpreting  model  results. 
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VEMTORA  COUNTY  FLOOD  ADVISORY 


PROVIDED  BY  THE 

CALXFORNIA-MEVAOA  RIVER  FORECAST  CENTER  OF  THE  NATIONAL  HEATHER  SERVICE 

r  ********************  A*******************  A  ******************  A*#*  A#*#**** 

FORECAST  PEAK  FLOWS  IN  THOUSAND  CFS  RESULTING  FROM  3  HOUR  PRECIPITATION 

3  HOUR  PRECIPITATION  (IN  INCHES) 


1 

2 

3 

4 

5 

SESPE  CREEK  NEAR  FILLMORE 

.45 

2.70 

10.89 

19.39 

27.98 

SANTA  PAULA  CR 

.12 

1.59 

5.35 

9.04 

12.72 

CALLEGOAS  CREEK  AT  CAMARILLO 

.23 

2.16 

8.88 

29.19 

50.90 

REV.  SLOUGH  (CAL  CK  FARM) 

.27 

.85 

4.02 

8.36 

12.82 

FAGAN  CANY(»( 

.02 

.09 

.36 

.54 

.71 

SANTA  ANA  CREEK 

.00 

.06 

.13 

3.98 

6.81 

COYOTE  CREEK 

.01 

.10 

1.37 

5.70 

9.47 

KATILIJA  CR  100%  BURNED 

.18 

.48 

6.11 

14.50 

22.52 

MATILIJA  CR  IF  UNBURNED 

.05 

.11 

.27 

2.39 

11.09 

ARROYO  SIMI  NR  SIMI 

.26 

.63 

9.95 

35.18 

55.67 

(Taylor  and  Weikel  1991) 


Figure  3.  FOeSCASTS  RASKD  ON  PBE-SHSOUj'IED  DEFTBS  OF  BVIUKK  RAINFAU. 


HnXKJLOGIC  HGDELIN5  -  GQIERAL 

A  large  number  of  hydrologic  models  are  available  for  developing  rain¬ 
fall  hyetografiis ,  simulating  runoff  and  determining  inundated  areas  for  his¬ 
torical  and/or  hypothetical  storm  events.  Such  models  can  be  used  effectively 
in  "ttie  development  of  hydrologic  criteria  for  flood  warning.  Characteristics 
of  two  such  models,  and  associated  data  management  software,  are  described  in 
the  following  sections.  The  programs  are  well  documented  and  can  be  used  on 
lEM-compatible  microccmiputers  as  well  as  UNIX-based  workstations.  Either 
metric  or  English  units  may  be  used. 


HEC-1.  Flood  Hvdrograph  Package 

Capabilities  -  The  fundamental  capability  of  HEC-1  (Hydrologic  Engineer¬ 
ing  Center,  1990a)  is  to  develop  discharge  hydrographs  for  historical  or  hypo¬ 
thetical  storm  events  at  specific  locations  in  a  basin.  The  basin  can  be 
subdivided  into  any  number  of  subbasins,  and  modeling  elements  such  as  uncon¬ 
trolled  reservoirs  and  diversions  can  be  accommodated.  The  program  includes 
options  to : 


o  optimize  values  for  unit  hydrogra^^  and/or  loss  rate  parameters 
o  optimize  values  for  routing  parameters 
o  simulate  enow  pack/snow  melt 
o  simulate  dam  overtopping/breaching 

o  incorporate  alteniative  land  use/development  conditions,  and  nul- 
tlple  storm  events,  in  a  single  program  application 

A  variety  of  alternative  methods  are  available  for  simulating  precipitation, 
losses,  base  flow,  runoff  transformation  and  routing. 

Precipitation  Computations  -  The  spatial  averaging  of  precipitation  can 
be  performed  externally  to  HEC-1  aixi  input  for  direct  use.  Alternatively, 
precipitation  for  individual  recording  and  non-recording  gages  can  be  speci¬ 
fied,  £dong  with  associated  wei^ting  factors  for  each  subbasln.  Additional 
wel^ting  to  accomnodate  gauge  bias  (e.g.,  difference  in  normal  annual  pre- 
clpltatlcxi  for  a  gauge  vs.  a  subbasin)  can  be  employed. 

Losses  and  fflrwa  -  Losses  can  be  specified  in  terms  of  (1)  an  in¬ 
itial  loss  and  constant  loss  rate,  (2)  a  four  parameter  exponent!^  loss  func¬ 
tion  (unique  to  HEX>*1),  (3)  an  initial  loss  and  a  "curve  number"  based  on  land 
use  and  DS  Soil  Conservation  Service  (SCS)  soil  classifications  (Soil 
Conservation  Service  1972),  (4)  the  Hoxtan  method,  and  (5)  the  (Srmen  and  Ampt 
method.  Base  flow  is  specified  by  means  of  three  input  variables. 

Runoff  Transform  -  Precipitation  excess  can  be  transformed  to  direct 
runoff  with  a  unit  hydrograph  or  kinematic  wave  techniques.  Unit  hydrograph 
options  allow  a  unit  hydrogra^^  to  be  input  directly  or  to  be  expressed  in 
terms  of  Clark,  Snyder  or  SCS  parameters.  The  kinematic  wave  option  permits 
deletion  of  subbasin  runoff  with  elanents  representing  up  to  two  overland- 
flow  planes,  two  collector  channels  and  a  main  channel. 

Routing  -  Primary  routing  options  are  the  Hiskingum-Cunge,  Modified  Puls 
and  Ikiekingum  methods.  For  the  Ikiskingum-Cunge  and  Modified  Puls  methods,  a 
routing  reach  can  be  specified  in  terms  of  a  length,  elope,  three  Manning  n 
values  (for  a  main  channel  and  left  and  ri^t  overbanks),  and  a  cross  section 
defined  with  ei^t  pairs  of  coordinates. 


HEC-2.  Water  Surface  Profiles 

HEC-2  (Hydrologic  Engineering  Center,  1991)  is  intended  for  calculating 
water  surface  profiles  for  steady,  gradually  varied  flow  in  natural  or  man¬ 
made  diannels.  Both  subcrltlcal  and  siQ)ercritical  profiles  can  be  calculated. 
The  effects  of  various  obstructions  such  as  bridges,  culverts,  weirs,  and 
structures  in  the  flood  plain  may  be  simulated.  The  computational  procedure 
is  based  on  the  solution  of  the  one-dimensional  energy  equation  with  energy 
losses  due  to  friction  evaluated  with  Manning's  equation.  The  energy  and 
energy- loss  equations  are  solved  iteratively  between  each  pair  of  cross  sec¬ 
tions  with  the  "standard  step"  method  (Henders<xi,  1966). 
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Data._StQrage  System  (HEC-D6S^ 


Data  management  Is  a  significant  aspect  of  hydrologic  evaluations.  The 
D6S  software  (Hydrologic  Engineering  Ceater,  1990)  is  intended  for  efficient 
management  of  time  series  data  of  any  l^pe,  sucdi  as  rainfall  hyetographs  and 
discharge  or  stage  hydrographs.  Falred-functlon  data,  such  as  stage-dlsdiarge 
rating  curves,  dlsduage-frequency  or  stage-frequency  relationships,  can  also 
be  aooonmodated.  The  D6S  system  Includes  a  set  of  utility  programs  Intended 
for  data  entry,  data  editing  and  graphic  displays.  Figure  4  Illustrates  the 
role  of  D6S  In  data  management.  A  typical  application  with  HBC-1  is  to  first 
store  observed  precipitation  and  dls<±arse  data  in  a  DBS  file  with  a  data 
entry  utility  program,  and  then  to  automatically  retrieve  such  data  as  part  of 
an  HEC-1  execution.  Simulation  results  can  be  written  to  the  same  D6S  file, 
and  another  utility  program  can  be  used  to  develop  grcphs  or  tabulations  of 
any  data  in  the  file. 


Figure  4.  BOLE  OF  DBS  IH  DATA  MANAnRWRNT 
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BBAL-TDIB  HYEBQLOGIC  MGDELINS 

SlmulatlcHi  iDodels  are  of  two  types,  eveat  and  contdnuouB.  Continuous 
modeling  generally  attempts  a  oontlnuous  accounting  of  soil  moisture,  whereas 
event  models  require  specification  of  initial  conditions  (e.g. ,  loss  rates) 
that  pertain  to  the  “event".  Advantages  of  event-type  approaches  are  that 
they  generally  do  not  require  representation  of  evapotranspiration  and  subsur¬ 
face  water  balances.  Because  event-type  approaches  are  simpler,  they  geansr- 
ally  use  fewer  paramefbers  and  are  easier  to  calibrate.  However  for  short-term 
forecasts,  there  may  be  substantial  uncertainty  with  respect  to  initial  condi¬ 
tions,  especially  after  dry  periods.  Tlie  following  sections  describe  first  an 
event- type  model,  HECIF,  and  then  a  continuous- type  model,  the  Sacramento 
Watershed  Model.  Comiimite  on  model  updating  are  also  provided. 


Event-type  Modeling 

Ihe  event-type  approadti  involves  the  following  steps  in  determining 
runoff  from  a  basin: 

o  specification  of  precipitation  (spatial  and  temporal  distribution) 
o  specification  of  “losses"  and  rainfall  excess 
o  transformation  of  rainfall  excess  to  direct  runoff 
o  specification  of  base  flow 

o  combining  of  base  flow  and  direct  runoff  to  obtain  total  runoff 

If  a  basin  is  divided  into  subbasins,  the  above  steps  are  performed  for  eacdi 
subbasin,  and  routing  and  combining  of  hydrographs  are  performed  as  required. 

HECIF,  which  is  an  adaptation  of  computer  program  HEC-1,  is  an  example 
of  an  event- type  model  used  for  forecasting  (Peters  and  Ely,  1985).  The  basic 
HEC-1  capabilities  for  calculating  runoff  with  a  unit  hydroep:aph  approach  from 
a  multi-subbasin  watershed,  and  for  parameter  optimizaticxi,  are  retained  in 
HECIF.  However,  HECIF  ocxitains  additional  capabilities  that  facilitate  the 
task  of  runoff  forecasting.  Aspects  of  application  of  HECIF  are  as  follows: 

1.  Forecasting  with  HECIF  is  intended  to  involve  a  "hands-on"  process  by 
which  the  analyst  can  readily  compare  simulated  hydrographs  with  ob¬ 
served  hydrogra^hs  (i;p  to  the  tiioe-of-forecast)  and  adjust  loss  rates, 
or  perhaps  other  parameters,  to  improve  results. 

2.  Forecasting  is  generally  performed  in  two  separate  executicxis  of  HECIF. 
In  the  first,  unit  hydrograph,  loss  rate  and  base  flow  parameters  are 
optimized  for  gauged  headwater  subbasins.  The  time  window  "T"  in 
Figure  5  is  the  period  over  vdiich  an  objective  function  to  optimize  the 
above  parameters  is  evaluated.  The  window  is  approximately  equal  to  the 
time  base  of  the  unit  hydrograph  for  the  subbasin.  An  objective  func- 
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a.  BASIN  b.  UNIT  HYOROQRAPH 


Figure  5.  n^RAHERBR  ESTIMmON  HTIH  HBCIF 


time  base  of  the  unit  hydrograph  for  the  subbasin.  An  objective  func¬ 
tion  is  minimized  by  a  univariate  gradient  technique  (Ford  et  al,  1980). 
The  objective  Sanction  is  as  follows: 


STDER  = 

objective  function 

ordinate  i  of  the  observed  hydrograj^i 
ordinate  i  of  the  ccmputed  hydrograph 
weighting  factor  applied  at  ordinate  i 
total  number  of  hydrograph  ordinates  encoiipassed 
by  the  objective  function 

Tte  equation  defining  the  wei^ting  factor  is  as  follows: 


there  STDER  = 

QCBSi  = 
QOOMPi  = 
WTi 

N  = 


= 

I 
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vrtiere  J  =  number  of  At  intervals  from  the  beginning  of  the  time 

period  for  parameter  estimation  (T)  to  the  time  of 
ordinate  i 

The  objective  function  is  a  quantitative  measure  of  the  goodness  of  fit 
of  the  calculated  hydrograph  to  the  observed  hydrograph.  Tlie  weighting 
factor  has  a  value  of  1  at  the  time-of -forecast,  and  to  a 

value  of  0  at  the  beginning  of  the  time  window  "T".  The  purpose  of  the 
weij^ting  is  to  Insure  a  relatively  close  fit  of  the  calculated  to  the 
observed  hydrograph  in  the  vicinity  of  the  time-of -forecast. 

The  optimization  process  has  built-in  constraints  that  prevent  physi¬ 
cally  unreasonable  values  for  the  parameters  to  be  optimized  (Hydrologic 
Engineering  Center,  1989).  For  example,  if  the  rainfall  is  oonoentrated 
very  near  the  tlme-of -forecast,  there  will  be  Uttle  hydrograph  “rise" 
with  tdilch  to  optimize  parameters.  In  this  case,  the  optimization  is 
permitted  only  for  base  flow  parameters. 

3.  Following  the  parameter  optimization  application  of  HBCIF,  the  analyst 
reviews  optimization  results  and  parameter  estimates  as  an  aid  to  set¬ 
ting  values  of  loss  rate  axxl  base  flow  parameters  for  the  remainder  of 
the  basin. 

4.  The  secxxid  application  of  HECIF  performs  runoff  computations,  and  rout¬ 
ing  and  combining  operations  throu^iout  the  basin.  At  each  location  for 
whldh  an  observed  hydrograph  is  available,  "blending"  can  be  perfoxmed. 

A  blended  hydrograph  consists  of  the  observed  hydrograph  vp  to  the  time- 
of-forecast  and  an  adjusted  simulated  hydrograph  eifter  the  tijne-of-fore- 
cast.  The  adjustment  is  made  either  by  a  vertical  shifting  of  the  simu¬ 
lated  hydrograph  with  a  constant  Increment  of  discharge  (positive  or 
negative),  or  by  providing  a  smooth  transition  from  the  observed  to  the 
unadjusted  simulated  hydrograph  over  six  time  intervals  following  the 
time-of- forecast .  The  transition  is  computed  by  linearly  diminishing 
the  "error"  (difference  between  the  observed  and  computed  discharge)  at 
ti,e  time-of- forecast  to  zero  over  the  six  time  intervals.  The  two  types 
of  blending  are  illustrated  in  Figure  6.  The  blended  hydrograph  is  used 
in  subsequent  routing  computations. 


HECIF  is  generally  used  in  oonjunctiai  with  the  program  PRECIP 
(Hydrologic  Engineering  Center,  1989)  for  processing  gauged  rainfall  data,  and 
DSS.  The  PRECIP  program  develops  hyetographs  of  spatial-average  rainfall  for 
each  subbasin  using  an  inverse  distance-squared  wei^tiiig  procedure.  If  data 
for  a  gauge  is  missing,  the  program  automatlcedly  obtains  data  for  the  next 
nearest  gauge.  The  DSS  graphics  utility  program  facilitates  data  review  and 
analysis,  as  well  as  evaluation  of  forecast  results.  Runoff  forecasts  based 
cxi  scenarios  of  future  (forecasted)  rainfall  can  be  readily  evaluated. 
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Figure  6.  THO  MEigCDG  FOR  EUHDIMB 


Ckaitinuous-tvc^  Modeling 

An  example  of  a  ODntinuous-type  model  used  for  real-time  applicatione  is 
the  Sacramento  Watershed  ttodel  (Bumash  et  al,  1973),  developed  by  the  Cali- 
fomia-Neyada  River  Forecast  Center  of  the  US  National  Weather  Service.  It 
has  been  in  use  for  a  number  of  years  by  that  agency  and  is  also  a  component 
of  the  National  Weather  Service  River  Forecasting  Syst^  (NWSRFS). 

The  Sacramento  Model  simulates  runoff  processes  in  headwater  basins.  It 
provides  a  conceptual  representation  of  (1)  soil  moisture  storage  (as  both 
tensicm  and  free  water)  at  two  levels,  an  upper  zone  and  a  lower  zone; 

(2)  direct  runoff  from  impervious  surfaces,  water  bodies  and  saturated  ground; 

(3)  percolation  from  the  upper  zone  free  water  storage  to  the  lower  zone; 

(4)  evaporation  from  surface  water  and  evapotranspi ration  from  tension  water 
storage;  (5)  interflow  from  upper  zone  free  water  storage;  and  (6)  baseflow 
and  subsurface  outflow  (out  of  the  basin)  from  lower  zcne  free  water  storage. 

A  unit  hydrograph  can  applied  to  surface  runoff  and  interflow,  and  the 
resulting  flow  can  optionally  be  routed  with  a  non-linear  "layered"  Muskingum 
method.  Figure  7  illustrates  storage  and  runoff  components  of  the  Sacramento 
Model. 


Figure  8  is  a  schematic  representation  of  analytical  aspects  of  the 
Sacramento  M^l.  Processes  are  characterized  in  terms  of  storages  of  speci¬ 
fied  capacities  v^ich  are  filled  from  precipitation  and  percolation,  and 
depleted  by  percolaticn,  evapotranspiration  and  lateral  drainage.  Table  1 
contains  brief  definitions  of  17  parameters  for  defining  the  production  of 
runoff  from  rainfall.  The  Table  does  not  include  parameters  associated  with  a 
unit  hydrograf*!  or  channel  routing. 
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(BumaBh  et  al  1973) 


Figure  8.  OOMPONOnS  OF  THE  SACRAHHtnX)  MODEL 
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Table  1.  BaraneterB  fbr  the  Sacramento  Nat^nsbed  Model 

PCTIM  Permanentiv  imoervioue  fraction  of  basin  continuous  wfth  stream  channels.  Rainfall  on  this  portion  of 
the  basin  bypasses  storages  and  contributes  instantiy  to  rurroff. 

AQIMP  Fraction  of  the  basin  which  becomes  impervious  when  upper-zorte  tension  storage  becomes  filled  (i.e., 
behavior  is  same  as  for  PCTIM). 

SARVA  Fraction  of  basin  covered  by  streams,  lakes,  and  riparian  vegetation  (generally  about  40%  to  100%  of 
PCTIM).  Evapotranspirabon  occurs  at  the  potential  rate  from  this  |M>rtion  of  the  basin. 

UZTWM  Maximum  depth  (over  norvimoervious  areas)  of  upper  zone  tension  water  storage.  This  storage  rnust 
be  filled  before  any  water  becomes  available  for  free  water  storage.  Water  from  this  storage  is  Tost*  by 
direct  evaporation  from  the  soil  surface  or  by  evapotranspiration  by  shallow-rooted  vegetation. 

UZFWM  Maximum  depth  (over  norv-impervious  areas)  of  u^/per  zone  free  water.  Water  from  this  zorre  feeds  the 
upper  zone  tension  water  storage,  percolates  to  the  lower  zone  or  corttributes  to  interflwv.  When  the 
free  water  storage  is  full,  rainfall  contributes  to  surface  runoff. 

UZK  Depletion  coefficient  for  upper  zone  free  Aoter  storage.  Upper  zone  lateral  drainage  (i.e.  interflow)  is  de¬ 
termined  on  a  daily  basis  as  the  product  of  this  coefficient  and  available  contents.  Percolation  from  the 
upper  zone  (during  a  time  irtferval)  is  acconrtmodated  prior  to  determination  of  interflow. 

PFREE  Proportion  of  water  percolating  from  the  iflsper  zone  to  the  lower  zone  which  passes  directly  to  the  free 

water  storages  without  contributing  to  the  tension  water  storage. 

LZTWM  Maximum  capacity  of  lower  zone  tension  storage.  Water  stored  here  can  only  be  removed  by 
evapotranspiration. 

L2FSM  Maximum  capacity  of  lower  zone  supplemental  free  water  storage.  Water  stored  here  contribittes  to 

baseflow  at  the  'supplemental*  rate. 

LZSK  Depletion  coefficient  for  lower  zone  supplemental  free  water  storage .  Supplemental  baseflow  (daily 
rate)  is  determined  as  the  product  of  this  coefficient  and  available  contents. 

LZFPM  Maximum  capacity  of  lower  zone  primary  free  water  storage.  Water  stored  here  contributes  to  baseflow 
at  the  ■primary*  rate. 

LZPK  Depletion  coefficient  for  lower  zone  primary  free  water  storage.  Primary  baseflow  (daily*  rate)  is  deter- 
rriined  as  the  product  of  this  coefficient  and  available  contents. 

RSERV  Decimal  fraction  of  lower  zone  free  water  storage  which  is  not  available  for  supplying  lower  zone  tension 
water  storage.  This  portion  of  the  lower  zone  free  water  storage  is  considered  to  be  below  the  root  sys¬ 
tem  that  reaches  to  tlie  lower  zone. 

SIDE  Decimal  fraction  of  baseflow  (entering  the  river)  that  exits  the  basin  without  passing  through  the  river 

channel.  For  example,  a  value  of  Q.5  indicates  that  a  quantity  of  baseflow  equal  to  50%  of  that  which  en¬ 
ters  the  river  leaves  the  basin  without  affecting  baseflow  in  the  river. 

SSOUT  Constant  streamflow  loss  term  representing  flow  through  porous  material  below  the  channel  bottom . 
This  term  is  commonly'  zero.  Rivers  in  glaciated  areas  may  require  use  of  this  parameter. 

2PERC  Proportional  increase  in  percolation  from  saturated  to  dry  condition. 

RBfP  Exponent  in  percolation  equation  which  determines  rate  at  which  percolation  demand  changes  from  the 
dry  condition,  (3^ERC  +  1)'PBASE,  to  the  wet  condition,  PBASE. 
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The  parameters  in  Table  1  provide  substantial  flexibility  for  repre¬ 
senting  rainfall- runoff  processes.  However  calibration  of  these  p6u:ametere 
can  be  a  very  challenging  task,  and  there  is  the  potential  for  mis-calibration 
such  that  the  model  is  treated  essentially  as  a  "black  box"  without  due  regard 
for  representation  of  the  essential  runoff  characteristics  of  the  basin.  A 
poorly  calibrated  model  is  ill-suited  for  prediction,  especially  where  condi¬ 
tions  differ  significantly  from  those  used  for  calibraticn.  It  is  generally 
necessary  to  have,  as  a  minimum,  several  years  of  continuous  precipitation  and 
streamflow  data  as  a  basis  for  calibration. 

A  simple  water  balance  of  annual  quantities  of  runoff,  rainfall  and 
evapotranspiration  over- a  period  of  several  years  can  be  used  to  gain  insist 
into  the  runoff  characteristics  of  a  basin.  Sudi  an  evaluation  can  facilitate 
recognition  of  subsurface  outflows  that  bypass  the  river  channel,  or  problems 
in  definiticxi  of  rainfall  volumes.  Runoff  (streamflow)  measurements  are 
generally  the  most  accurate  data  source  for  this  evaluation.  Spatially-aver¬ 
aged  rainfall  data  may  be  hi^ily  uncertain  because  of  large  sampling  errors 
(due  to  sparse  gage  networks)  and  because  of  errors  inherent  in  the  meas- 
uranent  of  rainfall.  Wind  effects  (e.g. ,  turbulence  over  the  gage  opening) 
tend  to  cause  rainfall  measuronents  to  be  biased  on  the  low  side,  as  illus¬ 
trated  in  Figure  9.  Underestimates  of  point  rainfall  depths  of  15%  or  more 
are  typical. 


(Bumash,  undated) 


Figure  9.  RAINFAU.  CATCH  VS.  HIND  SFEED 


14 


Evapotranspiration  must  be  estimated,  as  it  is  not  measured.  Actual 
evapotranspiration  can  differ  si®iificantly  from  evapotranspiration  potential, 
depending  on  water  availability.  Evapotranspiration  potential  is  scaoetimes 
based  on  application  of  monthly  or  seasonal  coefficients  applied  to  measured 
pan  evaporation.  However  there  is  much  uncertainty  in  this  approach,  as  the 
physical  processes  and  energy  fluxes  associated  with  pan  evaporation  can  be 
substantially  different  from  those  associated  with  evapotranspiration.  The 
Sacramento  Model  lumps  evaporation  from  water  bodies  and  moist  soil  with  eva¬ 
potranspiration  from  Vegetation.  Generally  the  latter  is  the  dominant  proc¬ 
ess.  Because  of  the  difficulty  in  using  pan  evaporation  data  as  a  surrogate 
for  the  total  baisin  evapotranspiration,  applications  of  the  Sacramento  Model 
commonly  utilize  direct  estimates  of  basin  evapotranspiration,  for  example 
using  data  such  as  that  provided  in  Figure  10  (Bumash,  undated). 
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Figure  10.  AVERAGE  DAILY  EVAPOTRANSPIRATION  DQIAND 
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Bumash  (Bumash,  undated)  suggests  initial  (typical)  values  for  the  17 
parameters  listed  previously.  He  also  provides  a  rationale  for  estimating 
values  for  several  of  the  parameters  from  direct  analysis  of  carefully  se¬ 
lected  portions  of  historical  rainfall  and  streamflow  data.  For  example, 
FCTIM,  the  Impervious  fraction  of  the  basin  contiguous  with  stream  channels, 
can  be  estimated  on  the  basis  of  small  runoff  events  following  Icaig  dry  peri¬ 
ods.  Presumably  surface  runoff  from  the  events  is  from  the  impervious  frac¬ 
tion  of  the  basin.  PCTIM  can  be  estimated  from  the  ratio  of  surface  runoff  to 
rainfall  accumulated  for  several  such  events. 

Because  of  interdependence  of  parameters  and  uncertainty  associated  with 
rainfall  data  and  evapotranspiration  estimates,  caution  must  be  exercised  in 
attempts  at  automated  calibration.  Bumash  recommends  review  of  monthly  error 
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sumnaries  with  due  cxOTsideration  of  uncertaijities  in  the  driving  inputs. 
Application  of  a  model  like  the  Sacramento  Model  requires  substantial  knowl¬ 
edge  of  the  runoff  characteristics  of  a  basin  and  skill  in  interpreting  cause 
and  effect  relationships  among  essential  processes.  Experience  in  applying  a 
model  of  this  type  can  itself  be  a  means  for  acquiring  understanding  of  basin 
rainfall-runoff  behavior. 


Model  Updating 

A  real-time  gauging  network  provides  observed  streamflow  data  up  to  the 
time-of-forecast.  This  data  can  used  to  adjust  model  inputs,  states  (i.e. 
volumes  of  water  in  storage)  and/or  parameter  values  so  that  the  model  is  more 
"in  tune"  with  current  conditions.  One  such  approach  is  to  apply  Kalman  fil¬ 
tering  for  automated  adjustment  (Georgakakos,  1986).  For  flash  floods,  auto¬ 
mated  adjustment  may  not  be  practical  because  by  the  time  a  stream  rise  oc¬ 
curs,  the  time  at  lAich  the  forecast  is  required  may  already  have  passed.  For 
an  event- type  model  like  HECIF,  model  updating  can  be  achieved  using  a  combi¬ 
nation  of  optimizaticoi  for  headwater  subbasins  and  manual  adjustment. 


MODEL  SEUBCmON  AMD  OSE 

Factors  to  consider  in  selecting  a  model  are  (1)  the  knowledge  and 
skills  of  the  model  user,  (2)  the  hydrologic  regime  to  which  the  model  will  be 
applied,  and  (3)  data  availability  and  data  requirements.  A  widely  cited 
study  by  the  World  Meteorological  Organization  (World  Meteorological 
Organization,  1975b)  indicates  that  in  humid  environments,  soil  moisture 
accounting  m^els  may  not  offer  much  advantage  over  simpler  "event"  approaches 
because  of  the  relatively  stable  moisture  conditions.  However  for  arid  and 
semi-arid  climates,  use  of  explicit  soil  moisture  accounting  models  (like  the 
Sacramento  Watershed  Model)  can  be  advantageous. 

Because  of  the  complex  nature  of  i^ysical  processes  that  constitute  the 
hydrologic  cycle,  the  heterogeneous  characteristics  of  a  watershed  (including 
the  subsurface),  and  the  uncertainty  associated  with  model  inputs,  rainfall- 
runoff  modeling  is  a  difficult  and  challenging  task.  The  skills  of  the 
analyst  applying  the  model  can  be  significantly  more  Important  than  the  model 
itself.  This  is  suppoi*ted  by  League  and  Freeze,  who  in  ccaiclusion  to  a  paper 
describing  a  comparison  of  modeling  techniques  or:  small  upland  catchments, 
state  the  following: 

"In  many  ways,  hydrologic  modeling  is  more  an  art  than  a  science, 
and  it  is  likely  to  remain  so.  Predictive  hydrologic  modeling  is 
normally  carried  out  on  a  given  catchment  using  a  specific  model 
under  the  supervision  of  an  individual  hydrologist.  The  useful¬ 
ness  of  the  results  depends  in  large  measure  on  the  talents  and 
experience  of  the  hydrologist  and  his  understanding  of  the  mathe¬ 
matical  nuances  of  his  particular  model  and  the  hydrologic  nuances 
of  his  particular  catchment.  ..."  (League  and  Freeze,  1985). 
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The  above  oonments  should  be  borne  in  mind  vdien  considering  the 
role  of  rainfall-runoff  modeling  in  flood  forecasting.  Modeling  can 
provide  valuable  information  vtiich,  Mhen  considered  in  relation  to  the 
currait  state  of  a  basin  and  meteorological  forecasts,  can  aid  in  the 
making  of  reasonable  flood  warning  decisions.  However,  there  will  al¬ 
ways  be  significant  uncertainty  associated  with  model  predictions,  and 
careful  interpretation  of  model  results  is  essential. 
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